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SIMPLIFIED SOLUTION FOR POINT CONTACT DEFORMATION 


BETWEEN TWO ELASTIC SOLIDS 

by David E. Brewe and Bernard J. Hamrock 

Lewis Research Center and 
U.S. Army Air Mobility R&D Laboratory 


SUMMARY 

A linear regression by the method of least squares is made on the geometric vari- 
ables that occur in the equation for point contact deformation. The ellipticity and tiie 
complete elliptic integrals of the first and second kind are expressed as a function of 
the X, y- plane principal radii. The ellipticity was varied from 1 (circular contact) to 
10 (a configuration approaching line contact) . The procedure for solving for these var- 
iables without the use of charts or a higjh-speed computer would be quite tedious. 

These simplified equations enable one to calculate easily the point- contact deformation 
to within 3 percent accuracy without resorting to charts or numerical methods. 


INTRODUCTION 

The classical Hertz solution requires the calculation of k, the ellipticity, and ^ 
and £, the complete elliptic integrals of the first and second kind, respectively. Pre- 
viously this entailed finding a solution to a transcendental equation that related k, ^ 
and £ to the geometry (ref. 1) . This was usually accomplished by some iterative nu- 
merical procedure (ref. 2) or with the aid of charts (refs. 1 and 3). 

In this report a linear regression by the method of least squares is made on these 
variables. The ellipticity and the complete elliptic integrals of the first and second 
kind are expressed as functions of the x, y-plane principal radii. The ellipticity was 
varied from 1 (circular contact) to 10 (a configuration approaching line contact). The 
resulting simplified equations enable one to easily calculate the point- contact deforma- 
tion to within 3 percent accuracy without resorting to charts or numerical methods. 



SYMBOLS 


A 

B 

E 


E’ 


£ 

£ 


F 

1 


R- 


X 


"y 

^Ax’ ’^Ay 


^Bx> ^By 



solid A 
solid B 

modulus of elasticity 



complete elliptic integral of the first kind 
£ expressed by method of least squares 
applied load 

complete elliptic integral of the second kind 

^ expressed by method of least squares 
ellipticity (ratio of semimajor to semiminor axis) 

k expressed by method of least squares 
effective radius of curvature in the principal x-plane 

effective radius of curvature in the principal y-plane 
principal radii of solid A 
principal radii of solid B 
curvature difference 

point- contact deformation at the center of contact 
point- contact deformation calciilated using curve fit data 
(l/Rjj + curvature sum 


POINT CONTACT DEFORMATION 

A widely used method of describing the geometry of two ellipsoidal solids, A and 
B, in contact (fig. 1) is to express it in terms of the effective curvatures, that is. 
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( 1 ) 
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R 



where 



(2) 


•'y "'Ay --By 

The variables and Ry represent the effective radius of curvature in the principal 
X and y planes, respectively. From reference 1 an auxiliary equation relating the 
curvature difference and the elliptic integrals can be written as 


^ / 2 f - .^(1 + D 

V ^(1-F) 


where 



( 4 ) 


( 5 ) 


Since ^ and are functions only of k, equation (4) indicates that the ellipticity k 

can be expressed strictly in terms of the principal plane curvatures. 

A numerical one-point iteration procedure (ref. 4) was used to generate k^ for a 

given (r_R„) for i = 1, 2,. . ., 39. The corresponding elliptic integrals were 
AX y/. 


evaluated munerically using a Landen transformation described by Bulirsh (ref. 5). 
Reference 6 indicates that k has been approximated as k = ^R Thus for a 


given set of pairs of data, 



2 .- • 



a power fit using a 


linear regression by the method of least squares resulted in the following equation; 


k = 


1.0339 



( 6 ) 
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The as 3 miptotic behavior of £ and ^ (ref. 7) was suggestive of the type of functional 
dependence that £ and migjit follow. As a resvilt, an inverse and logarithmic 
curve fit was tried for £ and respectively. The following expressions provided 
excellent curve fits: 


£ = 1 . 0003 + 


0. 5968 


(7) 


and 


^= 1 . 5277 + 0. 6023 In 

I'x 


( 8 ) 


Values of k, jf, and are presented in table I and compared with the numerically 
determined values of k, £, and All of the values of table I (except 

Ry/R^ = 1) were included in the least squares fit. An indication that the curve fits 
were excellent was that the coefficient of determination was not less than 0. 9997 for 
any of these expressions. A value of 1. 000 would indicate a perfect fit, while a value 
of zero would be the worst possible fit. 

From reference 1 the point contact deformation at the center of contact is 


where 



(9) 


( 10 ) 


If solid A is of the same material as solid B, then equations (9) and (10) can be com- 
bined to give 


2 3 

27T^E' V^R/ 


( 11 ) 


The point contact deformation can now be calculated in terms of the curve fit equations. 
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( 12 ) 


6 „ = 




Consequently, we see that for a given load and modulus, the deformation can be ex- 
pressed entirely in terms of and R^, Equation (11) can be determined also by 
using the charts by Jones (ref. 3). The methods of this report provide an alternative, 
simplified solution for the point- contact deformation, whereas Jones' charts are not 
always convenient or available. 

Table I compares 6p (determined from the curve fit equation) with 6p (deter- 
mined numerically). The ellipticity R was varied from 1 (a ball on a plate) to 10 (a 
configuration approaching line contact). For convenience the regression coefficients 
for k, £ , and ^ were rounded off to two places to the ri^t of the decimal. Conse- 
quently, a slight roimd-off error may add to or subtract from the percent error (per- 
cent given in table I). The percent error was determined as follows: 

(x - x)ioo 

6 =^^ — ^ ( 1 ^ 

X 


where 

X= Jk,jf, J^6pj 

X=-[k,J,?,6p} 

Note that for the case of a ball on a plate, which results in a circular contact, R^/R^ = 
1.0 and k=1.00. Further, for this case the elliptic integrals £ and reduce to 
7t/2. Consequently, 6p is determined directly without the need to determine k, £, 
and ^ via curve fitting, charts, or the numerical methods described in references 5 
and 6. The largest error (-5 percent) in determining 6p was for the circular contact 
case. If one eliminates the circular contact case, the largest error in determining the 
contact deformation using the curve fit values becomes less than 3 percent. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U. S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, May 18, 1976, 

505-04. 
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TABLE I. - COMPARISON OF THE NUMERICALLY DETERMINED VALUES WITH THE CURVE FIT 


VALUES FOR THE GEOMETRICALLY DEPENDENT VARIABLES 


Radius of 
curvature 
ratio, 

Ellipticity 

Complete elliptic 
integral of 
first kind 

Complete elliptic 
integral of 
second kind 

Point- contact deformation at center 
of contact, 6p 

k 

k 

J 

Percent 

error, 

e 

Determined 

numeri- 

cally, 

"p 

Determined 
from curve 
fit equation, 

6p 

Percent 

error, 

e 

£ 

£ 

Percent 

error, 

e 


sr 

Percent 

error, 

e 

^1.000 

1.00 

1.03 

3.00 

1.57 

1.60 


1.57 

1. 53 

-2. 55 

1. 230x10' 

1 . 168x10'^ 


2. 820 

1.99 

2.00 

. 50 

1.21 

1.21 

B 


2.15 

2. 15 

0 

1.020 

1.017 


5.314 

3.01 


-.33 

1. 11 

1. 11 

1 


2.53 

2. 53 

0 

.897 

.899 


8.330 

4.01 

4.00 

-.25 


1.07 

1 


2. 80 

2. 80 

0 

.814 

.816 

.25 

11.805 

4.99 

5.00 

.20 

1.05 

1. 05 

■ 


3.02 

3.01 

-.29 

.756 

.752 

-.53 

15.697 

5.97 

6.00 

. 50 

1.04 

1.04 

1 


3. 19 

3. 18 

-.25 

.706 

.701 

-.71 

19.971 

6.92 

7. 00 

1. 16 

1.03 

1.03 

■ 

■ 

3.33 

3.33 

0 

.667 

.662 

-.75 

24. 605 

7. 87 

8.00 

1.65 

1.02 

1.02 

■ 


3.46 

3.45 

-.24 

.636 

.628 

-1.26 

29. 576 

8. 80 

9.00 

2.27 

1.02 

1.02 

1 


3.57 

3.56 

-.22 

.608 

. 598 

-1.64 

34. 869 

9.72 

10.00 

2.88 

1.02 
1 

1.02 

1 


3.67 

3.66 

-.25 

.584 

. 571 

-2.23 


^Ry/Rj^ = 1. 000, then £ = tt/ 2 and k = 1. 00. There is no need to use the curve fit values for this case. 


/ 


F 



Figure 1. - Geometry of contacting elastic solids. 
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